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ABSTRACT 

The (i- and r-neutrinos emitted from a proto-neutron star are produced by nucleonic brems- 
strahlung NN — > NNvv and pair annihilation e + e~ — > z/p, reactions which freeze out at the 
"energy sphere." Before escaping from there to infinity the neutrinos diffuse through the "scatter- 
ing atmosphere," a layer where their main interaction is elastic scattering on nucleons vN — > Nv. 
If these collisions are taken to be iso-energetic as in all numerical supernova simulations, the neu- 
trino flux spectrum escaping to infinity depends only on the medium temperature Tes and the 
thermally averaged optical depth tes at the energy sphere. For tes = 10-50 one finds for the 
spectral flux temperature of the escaping neutrinos Xfl ux — 0.5-0.6 Tes- Including energy ex- 
change (nucleon recoil) in vN — > Nv can shift Tfl ux both up or down. AXfl ux depends on tes, 
on the scattering atmosphere's temperature profile, and on Tes- Based on a numerical study we 
find that for typical conditions ATfl ux /Tfl ux is between —10% and —20%, and even for extreme 
parameter choices does not exceed —30%. The exact value of ATfl ux /Tfl ux is surprisingly insensi- 
tive to the assumed value of the nucleon mass, i.e. the exact efficiency of energy transfer between 
neutrinos and nucleons is not important as long as it can occur at all. Therefore, calculating 
the Vp and v T spectra does not seem to require a precise knowledge of the nuclear medium's 
dynamical structure functions. 

Subject headings: diffusion — neutrinos — supernovae: general 



1. Introduction 

The role of [i- and r-neutrinos and anti- 
neutrinos in a supernova (SN) core is rather dif- 
ferent from that of v e and v e . After collapse, a 
huge amount of electron lepton number is trapped, 
leading to a large v e chemical potential, but no 
fi- or r- lepton number is present. The complete 
absence of r-leptons and the scarcity of muons im- 



plies that Vn 



and v T interact primarily by 



neutral-current processes while v e and v e also in- 
teract by more efficient charged-current reactions. 
The neutrino energies emitted from a SN core are 
a few tens of MeV, too low to produce muons 
or r-leptons, so that again v^ and v T interact 
only by neutral-current reactions with the matter 
above the SN core or in terrestrial detectors. On 
the other hand, v e and v e have charged-current 
reactions so that, for example, v e p — > ne + pro- 
vides the dominant SN neutrino signal in a water 



Cherenkov detector. Moreover, the v e and v e flux 
is thought to be responsible for re-heating the 
stalled shock in the delayed explosion scenario. 
Later, these neutrinos regulate the n/p ratio in 
the hot medium above the neutron star and thus 
govern the nucleosynthesis processes which take 
place in this region. Little wonder that in numer- 
ical studies far more attention has been paid to 
the treatment of v e and v e transport than to the 
other flavors. 

It is timely to take a fresh look at the v^ and 
v T spectra formation problem because the physi- 
cal interest in the flavor-dependent spectra is now 
much greater than it was a few years ago. The 
construction of a neutral-current neutrino obser- 
vatory is considered where the v^ and v T fluxes 
from a future galactic SN provide the dominant 
signal (Smith 1997; Boyd & Murphy 2001). Fur- 
ther, neutrino oscillations can partly swap the 
flavor-dependent spectra produced at the source, 
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especially if v e has a large mixing angle with the 
other flavors. A vast number of papers has been 
devoted to various consequences of flavor oscilla- 
tions on SN neutrinos and their detection (e.g. 
Raffclt f 996; Lunardini & Smirnov 2001), a topic 
of serious concern at a time when the phenomenon 
of flavor oscillations appears to be experimentally 
established. 

Moreover, there has been much progress in the 
numerical treatment of neutrino transport. New 
algorithms have been developed to effectively solve 
the Boltzmann collision equation in those critical 
SN regions where the neutrino mean free path is 
neither long nor short relative to the important ge- 
ometric scales (Burrows et al. 2000; Mezzacappa 
et al. 2001; Rampp & Janka 2000). With vastly 
increased CPU power one is beginning to perform 
simulations where a reliable transport scheme is 
self-consistently coupled with the hydrodynamic 
evolution so that the calculated multi-flavor neu- 
trino fluxes and spectra will depend only on the 
adopted input physics. 

Our present approach is complementary to 
these global numerical simulations. We study v^ 
and v r transport and spectra formation in the 
framework of the simplest possible model that 
incorporates enough of the essential physics to 
mimic the full problem. This approach allows us 
to ascertain the significance (or irrelevance) of mi- 
croscopic input-physics variations that may mod- 
ify the spectra. The nontrivial insights gathered 
from this study can serve as a basis for a more in- 
formed choice about the micro-physics that should 
be implemented in a full simulation. In addition, 
our approach has the pedagogical benefit of pro- 
viding a simple and transparent framework for 
understanding the crucial physics. 

The spectra formation problem is schematically 
illustrated in Fig. 1. The electron (anti-)neutrinos 
are kept in thermal equilibrium by beta processes 
up to a radius usually referred to as the "neutrino 
sphere." Beyond this radius the neutrinos stream 
off freely, their spectrum representing the medium 
temperature at the neutrino sphere. Of course, 
this picture is crude because the interaction cross 
section varies as e 2 (neutrino energy e) so that 
different energy groups decouple at different radii 
and thus at different medium temperatures. 

The other flavors interact with the medium pri- 
marily by neutral-current collisions on nuclcons 
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Fig. 1. — Schematic picture of neutrino spectra 
formation in the atmosphere of a SN core. 



vN <-> Nv, a reaction which is sub-dominant 
for the electron flavor. The nucleon mass m = 
938 MeV is much larger than the relevant temper- 
atures which are around T = 10 MeV so that en- 
ergy exchange between neutrinos and nuclcons is 
inefficient. However, nucleon-nucleon bremsstrah- 
lung NN <-> NNvv as well as the leptonic pro- 
cesses e + e~ <-> vv and ve <-> ev allow for the ex- 
change of energy and the creation or destruction 
of neutrino pairs and thus keep neutrinos in local 
thermal equilibrium up to a radius where these re- 
actions freeze out, the "energy sphere." However, 
the neutrinos are still trapped by vN <-> Nv up to 
the "transport sphere" whence they stream freely. 
Between the energy and transport spheres, neutri- 
nos propagate by diffusion. This region plays the 
role of a scattering atmosphere. 

In all numerical simulations of SN neutrino 
transport the neutrino collisions in the scatter- 
ing atmosphere were treated as iso-energetic so 
that the energy ei of the outgoing neutrino in 
vN — ► Nv was set equal to the energy t\ of the 
initial state. The main motivation for this approx- 
imation was its numerical simplicity and the lack 
of a compelling interest in details of the emerg- 
ing v^ and v T spectra. It is clear, however, that 
iso-energetic collisions are not a particularly good 
approximation. In Fig. 2 we show the distribution 
of final-state energies £2 when e\ — 30 MeV and 
the medium temperature is 10 MeV. A typical 
nucleon velocity is then about 20% of the speed 
of light so that it is not surprising that even after 
a single collision the neutrino energy is consider- 
ably smeared out. Since neutrinos interact many 
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Fig. 2. — Distribution of final-state energies €2 
of a neutrino with initial energy e\ = 30 MeV, 
scattering on non-degenerate nucleons in thermal 
equilibrium with T = 10 MeV. Details of how to 
calculate this plot are described in Appendix B.l. 



times in the scattering atmosphere, and since the 
medium temperature decreases between the en- 
ergy and transport spheres, there can be a signif- 
icant downward adjustment of the neutrino ener- 
gies (Janka et al. 1996; Hannestad & Raffelt 1998). 
The main purpose of the present paper is to pro- 
vide a conceptual understanding and a quantita- 
tive estimate of the magnitude of this effect. 

To address this problem we simplify the model 
of Fig. 1. The very concept of an "energy sphere" 
suggests that one should think of it as a source 
of thermal neutrinos which subsequently diffuse 
through the scattering atmosphere. Taking this 
concept literally amounts to the simple picture 
illustrated in Fig. 3. One no longer worries about 
detailed processes like NN bremsstrahlung to 
thermalize the neutrinos, but directly feeds a ther- 
mal flux into the scattering atmosphere. 

Section II of our paper is devoted to showing 
that this simple picture actually provides a sur- 
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Fig. 3. — Schematic picture of our simplified treat- 
ment of the scattering atmosphere. Tes is the 
medium temperature at the energy sphere. 



prisingly accurate representation of the spectra 
formation problem. The neutrinos streaming off 
the transport sphere then have fluxes and spectra 
which depend only on the temperature Tes and 
the thermally averaged transport optical depth 
tes a t the energy sphere which here coincides with 
the bottom of the scattering atmosphere. 

As a next step in Section III we study a scatter- 
ing atmosphere with a blackbody boundary con- 
dition at the bottom and with iso-energetic vN 
collisions as the only neutrino interaction channel. 
We derive an explicit relationship between Tes and 
the spectral flux temperature Tfl ux of the escaping 
neutrinos as a function of tes- Comparing with 
full-scale numerical simulations indicates that this 
exceedingly simple model accounts for the main 
features of the and v T spectra. 

Then in Section IV we include nucleon recoils 
in this model. We consider different types of tem- 
perature profiles to estimate the shift of the flux 
temperature and identify the critical parameters 
which govern ATfl ux . 

In Section V we summarize and discuss our 
findings. Many technical details, especially re- 
garding our implementation of neutrino-nucleon 
interactions with recoil energy transfer and NN 
bremsstrahlung, are documented in a series of Ap- 
pendices. 

2. Energy Sphere as a Blackbody Surface 

The neutrino sphere, transport sphere, or en- 
ergy sphere are only approximate concepts be- 
cause of the energy dependence of the neutrino 
cross sections. Still, we presently argue that the 
concept of the energy sphere as a well-defined 
blackbody surface at the bottom of the scattering 
atmosphere is much better suited to understand 
the neutrino spectra than one may have hoped. 

To get started we study a practical example. 
We have performed a and v T transport calcula- 
tion with the simple Monte Carlo code described 
in Appendix E. We use plane-parallel geometry 
with the medium temperature and density follow- 
ing the power laws 
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and T = T ( 



where for the present example p = 10 and q 
10/4. The other parameters are r = 10 km, p 
2 x 10 14 g cm" 3 and T = 31.66 MeV. 
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The neutrinos interact by iso-energetic vN scat- 
tering and by NN bremsstrahlung. This choice of 
micro-physics is orthogonal to most previous stud- 
ies which ignored bremsstrahlung, but included 
e + e~ «-> vv and ve <-» ev. In this traditional 
approach spectra formation is conceptually even 
more complicated because the number-changing 
process e + e~ <-> i/z/ freezes out more deeply than 
the energy-changing process ve <-> so that in 
addition to the transport and energy spheres there 
is a "number sphere" (Suzuki 1989). Later it 
was recognized that NN bremsstrahlung is more 
important than e + e~ annihilation (Suzuki 1991, 
1993). Hannestad & Raffelt (1998), Thompson, 
Burrows & Horvath (2000), and Burrows et al. 
(2000) confirm the importance of N N bremsstrah- 
lung. Moreover, Hannestad & Raffelt (1998) find 
that the freeze-out spheres for bremsstrahlung and 
for ve <-» ev roughly coincide in their example of 
a SN core. Burrows et al. (2000) find that the 
v^ spectrum depends sensitively on the assumed 
strength of NN bremsstrahlung, suggesting that 
v^ thcrmalization is not vastly dominated by ve 
scattering. 

Therefore, reducing the thcrmalization pro- 
cesses to NN bremsstrahlung likely captures the 
dominant effect. As this process allows for both 
the change of neutrino number and energy, the 
concept of an energy sphere or "thcrmalization 
depth" is unambiguous, but of course energy de- 
pendent. Ignoring ve scattering has the additional 
benefit that it becomes self-consistent to ignore 
the chemical composition of the medium. While 
vN scattering and NN bremsstrahlung do depend 
on the chemical composition, this effect is weak 
for vN scattering. While it may be strong for NN 
bremsstrahlung, the existing calculations of this 
process are uncertain to perhaps a factor of 2 (see 
Appendix C) so that one would not gain much by 
modeling the composition dependence. 

Once the micro physics has been fixed it is 
straightforward to locate the energy sphere by 
virtue of a classic argument. Assume the neutri- 
nos have only two interaction channels, i.e. iso- 
energetic scattering with the inverse transport 
mean free path (mfp) A^ 1 and a reaction with 
the mfp A^ 1 which changes the energy by a large 
amount and thus leads to quick thermalization. In 
this situation the optical depth J dr A^ 1 for energy 
exchange is not relevant because neutrinos with 



a large A^ are trapped and thus have a greater 
chance of participating in an energy-changing re- 
action. This reasoning leads to the optical depth 
for thermalization (Shapiro & Teukolsky 1983) 

TthermM = /"V ^V) V) + *J (*")] ■ 

J r 

(2) 

The radius tes of the energy sphere is implied by 

2 

Ttherm(rEs) = g- (3) 

Averaging r t horm over a thermal neutrino spec- 
trum at the local medium temperature leads in our 
example to tes = 13.14 km for the location of the 
energy sphere with a temperature Tes = 16 MeV. 

In order to discuss non-equilibrium neutrino 
fluxes and spectra we need to characterize the neu- 
trino distribution function by a few intuitive pa- 
rameters. One is the spectral temperature which 
we define as 

3 3 jo°° rfe S-i 

Here, /(e,/u) is the neutrino distribution function 
with e the energy and \i the cosine of the angle 
between the neutrino momentum and the radial 
direction. (Even for our plane-parallel geometry 
we will often say "radial" when we mean perpen- 
dicular to the medium layer.) We always ignore 
Pauli blocking, i.e. the Fcrmi-Dirac nature of neu- 
trinos. The Boltzmann collision equation is then 
linear in /, a considerable numerical simplification 
which is justified because the v^ and v T distribu- 
tions do not have chemical potentials in a SN core. 
Moreover, in the scattering atmosphere their dis- 
tribution function is diluted, further reducing the 
relevance of phase-space blocking. For our study 
the simplification of ignoring blocking effects far 
outweighs the small loss of precision. Boltzmann 
statistics implies that in equilibrium neutrinos fol- 
low the distribution function / = e~ e / T so that 
(e) = 3T. Inverting this relationship for the non- 
equilibrium case leads to our definition of T. t . 

In Fig. 4 we show the T* profile for our ex- 
ample as a solid line. (The slight wiggles arise 
from the limited statistics of the numerical Monte 
Carlo calculation.) We observe that T* coincides 
with the medium temperature T me( j up to the en- 
ergy sphere; for larger radii the two curves sepa- 
rate. T* still drops considerably until it reaches 
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an asymptotic value which is about half of the 
energy-sphere temperature Tes- This behavior is 
understood by the e 2 dependence of the scattering 
cross section which causes higher-energy neutrinos 
to be trapped more effectively than lower-energy 
ones. Hence, the escaping neutrino flux must be 
biased towards lower energies. 

Therefore, one should carefully distinguish be- 
tween the spectral temperature as defined in 
Eq. (4) and the spectral flux temperature 

r+i 



Tflux = 



_ Jo°° dee Li d(J,/J,f{e,fi) 



3 JIT de S-i /( e ' 



(5) 



This is (e)/3 of those neutrinos which actually 
stream. Note that apart from overall factors 
2tt d/j, /j, f(e, fj,) represents the particle flux. 
Tflux does not include the isotropic part of the dis- 
tribution which determines when the neutrinos 
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Fig. 4. — Numerical example for and v T trans- 
port as described in the text. The inner boundary 
condition was at 10 km. The thin dashed lines 
in the upper panel refer to the case without TV TV 
bremsstrahlung where a blackbody boundary con- 
dition with T ES = 16 MeV is used at the energy 
sphere. 



are trapped. At large distances where neutrinos 
stream freely we have f(e, fj) = for backward 
directions so that T* and Tfl ux are identical up 
to small differences caused by the exact angular 
distribution. 

In Fig. 4 we show the Tfl ux profile for our nu- 
merical example. At large radii indeed T* w Tfl ux . 
Moreover, Tfl ux stays constant throughout the 
scattering atmosphere because no energy exchange 
with the medium is possible. Flux conservation 
in the scattering atmosphere is apparent from the 
lower panel where we show the particle and energy 
fluxes which both remain constant from roughly 
the energy sphere outward. 

It appears that the separation of the T* and 
Tmed profiles, as well as the saturation of Tfl ux and 
the particle and energy fluxes all happen within a 
narrow range of radii around the energy sphere. 
This observation motivates us to repeat the simu- 
lation with the lower boundary now at tes > switch- 
ing off NN bremsstrahlung entirely, and injecting 
neutrinos with a thermal spectrum corresponding 
to Tes = 16 MeV, i.e. enforcing a neutrino black- 
body surface at r^s- The resulting temperature 
profiles are shown in Fig. 4 as dashed lines. Tfl ux 
is now strictly horizontal. Likewise, the energy 
and particle fluxes are now strictly conserved so 
that their profiles are horizontal lines not shown 
in this plot. 

The T» and Tfl ux profiles arc surprisingly close 
to the previous case where neutrinos thermalized 
by NN bremsstrahlung. Thus, the concept of an 
energy sphere as a blackbody surface looks encour- 
agingly well justified. 

One reason for the extraordinary "sharpness" 
of the energy sphere in our example is the un- 
usual behavior of the neutrino absorption rate by 
inverse NN bremsstrahlung; Eq. (C12) implies ap- 
proximately A^ r g ms oc e _1 . This is one of the few 
neutrino reaction rates which decrease with en- 
ergy; the scattering rate A^ 1 increases as e 2 . With 
+ A^ rems w A^ 1 we find that the effective ther- 
malization rate under the integral in Eq. (2) varies 
as e 1 / 2 . Therefore, different energy groups decou- 
ple thermally at similar radii. 

Another reason for a sharply defined energy- 
sphere temperature is that the density as a func- 
tion of radius falls much more steeply than the 
temperature. Therefore, as a function of optical 
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depth the temperature varies slowly. In our ex- 
ample A^ 1 oc p and approximately A^ rems oc p 2 T 3 

so that (A^Ay 1 ) 1 / 2 oc p 3 / 2 T 3 / 2 or r th crm oc 
( ro / r )3(p+<jT/2-i SQ that T a T » orm with w = 

2q/[3(p + q) — 2]. In our explicit example this 
is T oc t^J^. Therefore, the energy spheres for 
the different energy groups are at similar temper- 
atures. 

We thus expect that a blackbody boundary 
condition at the bottom of the scattering atmo- 
sphere provides us with a reasonably accurate un- 
derstanding of the i/fj, and v r spectra formation 
problem. On the other hand, for the problem 
of hydrostatically self-consistent atmospheres with 
multi-flavor neutrino transport (e.g. Schinder & 
Shapiro 1983) it would be physically more appro- 
priate to use a prescribed neutrino energy flux 
from the inner star; the atmosphere would then 
be allowed to adjust in order to transport exactly 
this flux. In our study, on the other hand, we al- 
ways assume a fixed background model where the 
neutrino fluxes are determined by externally pre- 
scribed density and temperature profiles as well 
as the physical assumptions about their interac- 
tion with the medium. In that sense our study 
is similar in spirit to the Monte-Carlo studies of 
Janka & Hillcbrandt (1989a, b) or the more recent 
work of Burrows et al. (2000). We do not wish to 
imply that in a real SN core the neutrino fluxes 
are determined by the temperature at the base of 
the scattering atmosphere unless this temperature 
corresponds to a self-consistent stellar structure. 

3. Spectral Modification by the Scattering 
Atmosphere 

3.1. Flux Transmission 

In our simple picture where the scattering at- 
mosphere is irradiated by a thermal neutrino flux 
from below, part of that flux will be reflected, part 
of it transmitted. The reflected flux will be re- 
absorbed because a blackbody surface is by defi- 
nition a perfect absorber. To determine the emerg- 
ing and v T flux it is therefore enough to calcu- 
late the energy-dependent flux transmission fac- 
tor, i.e. the fraction of the primary flux that is 
transmitted without being reflected back into the 
radiating surface. 

For sufficiently low energies the scattering at- 
mosphere is transparent so that the transmission 



factor is unity. In the opposite limit where the 
neutrino mfp is short compared with the geometric 
dimension of the layer (diffusion limit), the Boltz- 
mann collision equation can be solved analytically 
(Appendix A), leading for a plane-parallel geome- 
try to a transmission ratio of 



s(tes) 



3tes 



(6) 



Here, t es is the transport optical depth of the en- 
ergy sphere 



/ + 1 i 
dcosO (1 — cos 9) 
-i 



dcos9 
(7) 

where tib is the number density of scatterers 
("baryons"), da jd cos 9 the differential scattering 
cross section, and r^s the lower boundary of the 
scattering atmosphere, i.e. the energy sphere. 

In the intermediate regime tes = 0(1) we 
solve the Boltzmann collision equation numeri- 
cally. To do so we have to commit ourselves to 
a specific form of the differential scattering cross 
section. We use the "dipole formula" da jd cos 9 = 
a (l + 6cos6*)/2 where — 1 < b < +1 and the total 
cross section is a^. When expressed in terms of the 
transport optical depth tes, the flux transmission 
is independent of b for small and large tes- For 
intermediate tes the variation with b is less than 
±1%, i.e. for all practical intents and purposes it 
is enough to use the transport optical depth as the 
only relevant parameter. 

The numerical transmission factor for b = is 
shown in Fig. 5; for other b- values it is within the 
line-width of this curve. An analytic approxima- 
tion is 



s(tes) = 
x 1 



1 + 



3 T E s 



0.033 



1 + 1.5 (£ T + 0.17) 2 + 0.5 (£ T + 0.32) 6 

(8) 

where £ T = log 10 ('TEs)- The first bracket alone 
overestimates the numerical result by less than 4% 
in the entire range < tes < oo. The second 
bracket improves the approximation to better than 
±0.05%. 

For a spherically symmetric geometry these re- 
sults need to be modified. In Appendix A we show 
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Fig. 5. — Flux transmission factor for a plane- 
parallel scattering atmosphere with transport op- Fi S- 6 -— Transmitted flux for the indicated values 

tical depth tes • ofrES - 



that in the diffusion limit the transmission factor is 
the same as for a plane-parallel geometry, provided 
one includes a factor (tesA) 2 under the integral 
in Eq. (7). The resulting quantity r^g, of course, 
is no longer the transport optical depth, but for 
the purposes of flux transmission it plays an anal- 
ogous role. For the spherical case we have not 
investigated the intermediate regime Tg S = 0(1), 
but assume that the plane-parallel interpolation 
formula Eq. (8) applies with reasonable accuracy. 

3.2. Spectral Modification 

Next we study the flux and spectrum of the 
transmitted neutrinos if the primary flux is ther- 
mal at a temperature Tes- We characterize the 
"thickness" of the plane-parallel scattering at- 
mosphere by tes, defined as the transport op- 
tical depth at the energy sphere, averaged over 
a Maxwell-Boltzmann spectrum at Tes- Tak- 
ing the scattering cross section to vary as e 2 we 
note that (e 2 ) = 12T| S so that tes amounts to 
the transport optical depth at the fixed energy 
e = VV2T ES w 3.4641 T ES . Put another way, 
we use tes (e) = (f E s/12)(e/r E s) 2 to calculate the 
energy-dependent transmission factor. 

In Fig. 6 we show the energy-dependent trans- 
mitted flux for several values of tes • One parame- 
ter to characterize these spectra is their flux tem- 
perature Tfl ux as defined in Eq. (5). In Fig. 7 we 
show the same spectra after re-scaling the horizon- 
tal axis with the appropriate Tfl ux and normalizing 
them. This representation shows how closely these 
functions mimic true thermal spectra. 
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Fig. 7. — Transmitted flux (normalized) for the 
indicated values of tes where the neutrino energy 
is in units of Tfj ux . 



The flux temperature Tfl ux in units of Tes as 
a function of tes is shown in the upper panel 
of Fig. 8. For tes <C 1 we have, of course, 
Tflux = Tes- In the opposite limit tes — » oo 
the transmission factor is proportional to e~ 2 so 
that the transmitted spectrum is proportional to 
e -e/T E s ra ther than the primary e 2 e _e / TES , 
: ?es • Th 



imply- 
ing Tfl ux = i Tes ■ Tflux is significantly smaller 
than T E s even f° r relatively small values of f E s- 
An analytic approximation formula is 



Tflu 



Tes 



+ 



3 ' 3(1 + 0.35 fgjp) 
1 -0.03(1- 0.75 y 2 )e-°- 27 ^ 



(9) 
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where y — log 10 (fEs) — 0.9. The first bracket alone 
approximates the true result to better than 4% 
while the second bracket improves the approxima- 
tion to better than ±0.3% everywhere. To study 
power-law models it is also useful to have a power- 
law representation. We find that 



Tes 



0.83 

,==0.128 
T ES 



(10) 



works well in the range 1 < tes J$ 100 (dashed 
line in the upper panel of Fig. 8). 
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Fig. 8. — Spectral characteristics of the transmit- 
ted flux. 



The deviation of the spectrum from a thermal 
shape may be characterized by the "pinching pa- 
rameter" 

>-\% 

where p = 1 for a Maxwell-Boltzmann spectrum. 
If the high- and low-energy parts of the spectrum 
are relatively suppressed, the spectrum is called 
"pinched" (p < 1), in the opposite case "anti- 
pinched" (p > 1). The limiting cases are p = 1 
for tes - * and | for tes — » oo. Intermediate 
values are shown in Fig. 8 and Table 1. 

The concept of spectral pinching refers to the 
width of the energy distribution. Therefore, an- 
other choice for the pinching parameter might 
have been something like ((e 2 ) - (e) 2 )/(e) 2 . This 
definition is equivalent to p aside from a shift of 
the zero point and the normalization. Either way 
one uses the quadratic moment (e 2 ) in addition to 
(e) to characterize the spectrum. 

In the literature, non-thermal SN neutrino 
spectra are often described by a degeneracy pa- 
rameter r\. This representation has the disadvan- 
tage that anti-pinched spectra are not covered be- 
cause the limit i] — > — 00 represents a Boltzmann 
distribution. 

A further characteristic of the transmitted flux 
is its luminosity. It may be larger or smaller than 
given by the Stefan-Boltzmann law in terms of 
Tflux- We call the ratio of the true luminosity rela- 
tive to this Stefan-Boltzmann value the "greyness 
factor" (p. When <fi < 1 the transmitted flux is 
"grey" in that it carries less energy than a black- 
body flux of the same spectral shape. In the bot- 
tom panel of Fig. 8 we show </>; some values arc 
given in Table 1. 

If the scattering atmosphere is very thick 
(tes ^ 1) the emerging spectrum approaches a 
limiting spectral temperature of §7es as explained 
above Eq. (9), while the flux is arbitrarily diluted. 
Therefore, the greyness factor will become arbi- 
trarily small. It is remarkable, however, that 0=1 
within ±20% for tes as large as 30. Therefore, if 
in realistic examples of density and temperature 
profiles the energy sphere is at tes < 30, then 
the emerging flux will be surprisingly close to 
a Stefan-Boltzmann flux as calculated with the 
spectral temperature Tfl ux . Put another way, for 
tes < 30 the overall flux dilution and spectral 
shift conspire to mimic the Stefan-Boltzmann law 



8 



Table 1 

Spectral modification of neutrino flux. 



tes 


7flux/?ES 


V 


4> 





1 


1 


1 


1 


0.82 


0.99 


1.18 


3 


0.72 


1.01 


1.17 


10 


0.62 


1.05 


1.04 


30 


0.54 


1.10 


0.79 


100 


0.47 


1.17 


0.49 


oo 


1/3 


3/2 






for the emerging flux. This unintuitive behav- 
ior is a consequence of the e 2 dependence of the 
transport cross section. If the cross section were 
energy independent, then the flux would be di- 
luted without a spectral shift, and even for small 
tes the emerging flux would appear diluted (grey) 
relative to its spectral temperature. 

The neutrino fluxes and spectra from SN cores 
are usually described as showing equipartition of 
energy between different flavors, yet a hierarchy 
of spectral temperatures. At late times when the 
star is settled and the density and temperature 
profiles in the atmosphere are steep, the radiating 
surfaces for different flavors arc similar. There- 
fore, as stressed by Janka (1995), similar luminosi- 
ties with different spectral temperatures must be 
explained by the flux dilution caused by the scat- 
tering atmosphere. Janka (1995) did not worry 
about the spectral shift caused by the e 2 variation 
of the transport cross section so that his conclu- 
sions would have applied to any optical depth of 
the energy sphere. Our more refined treatment 
suggests that large spectral differences and similar 
luminosities for the different flavors at late times 
require average optical depths of the energy 
sphere exceeding about 30. However, it is precis- 
esly at late times when Y e is small and the den- 
sity gradients are steep that the NN bremsstrah- 
lung process is likely to dominate over the leptonic 
processes. Once NN bremsstrahlung is included 
in future numerical simulations it remains to be 
seen if at late times the energy sphere moves to 
such large optical depths as to cause significantly 
"grey" or v T fluxes. 



In summary, for values of tes up to 30 the 
transmitted flux resembles a thermal spectrum in 
both its absolute flux as well as its spectral shape. 
However, the relevant temperature is not Tes, but 
rather a reduced temperature Tfl ux . For tes = 10— 
30 we find T flux /T ES = 0.62-0.54. Therefore, we 
expect that the ty, and v T flux from a SN core 
is roughly thermal with a temperature Tfl ux w 
0.6 Tes and a luminosity given approximately by 
the Stefan-Boltzmann law, taking Tfl ux for the 
temperature and the area of the energy sphere 
for the radiating surface. Only for tes < 30 is 
the emerging flux significantly diluted below the 
Stefan-Boltzmann level corresponding to its spec- 
tral temperature Tjj ux . 

3.3. Comparing with Full-Scale Simula- 
tions 

3.3.1. Janka & Hillebrandt (1989b) Model 

As a "reality check" on our explanation of 
and v T spectra we turn to comparing our pre- 
dictions with full-scale numerical simulations in 
the literature. As a first case we consider the 
Monte Carlo transport simulations of Janka & 
Hillebrandt (1989a, b) where the multi-flavor neu- 
trino fluxes were calculated on the basis of a fixed 
background model. This model consisted of ten 
radial zones, each 12 km wide; in Fig. 9 we show 
the temperature and density profiles. (I thank 
Thomas Janka for providing these profiles.) The 
local and v T temperature T* found in this 
Monte Carlo simulation is indicated by filled cir- 
cles in the third panel of Fig. 9; the flux temper- 
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Fig. 9. — SN model for the Monte Carlo transport 
study of Janka & Hillebrandt (1989b), Model II. 
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Fig. 10. — SN model 324 ms after bounce from 
Newtonian calculation (Messer et al. 2001). 
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ature Tfl ux by open circles. 1 As expected, T* and 
Tflux coincide for large radii, but are vastly differ- 
ent in regions where the neutrinos still scatter. 

In the second panel of Fig. 9 we show the av- 
erage transport optical depth f . To this end we 
calculate the transport optical depth, modified for 
the spherical case, and average it thermally with 
the local medium temperature. For the transport 
opacity we use Eq. (B3) which ignores the vector- 
current contribution and also ignores nucleon de- 
generacy effects, two approximations which work 
in opposite directions, and anyway, are small cor- 
rections for the situation at hand. 

On the basis of this f(r) and T me <i(r) we cal- 
culate for each radius the predicted flux temper- 
ature Tfl ux (rEs) of the escaping neutrino flux if 
the energy sphere were at that radius. This is the 
lowest curve in the third panel of Fig. 9. The ac- 
tual Tflux at large radii found in the Monte Carlo 
study is shown as a horizontal dotted line in the 
third panel. It intersects the T nux (rEs) curve at 
a radius of about 50 km, indicated by a vertical 
dotted line. This radius is tes if our picture is 
correct. We note that at r^s — 50 km we have 
f w 10. 

The neutrino temperature T* (top curve in the 
third panel) separates from T moc i at about the 
same radius so that it is consistent to interpret 
50 km as ves- 

Likewise, the numerical particle flux shown in 
the fourth panel reaches its asymptotic value at 
about the same radius so that, indeed, the energy 
sphere must be around 50 km. Within our sim- 
ple picture we can also predict the particle flux, 
assuming that the energy sphere is at a given ra- 
dius. This prediction is shown as a thin line in 
the fourth panel. If tes = 50 km we would pre- 
dict that a neutrino flux of about 3.5 x 10 57 s _1 
is emitted from this SN model. Comparing this 
with the Monte Carlo result of 2.6 x 10 57 s _1 , the 
agreement is not too bad. 

In this numerical study the only number- 
changing reaction was e + e~ <-> vD, for energy 
exchange there was also ve <-» ev. As these pro- 



The data are extracted from Table IV of Janka & Hillc- 
brandt (1989b). T„ is obtained by dividing in the 

fifth column by 3.15 rather than by 3 because Fermi-Dirac 
statistics were used here. Tfl ux is obtained by multiply- 
ing T* with £ e /£v which are given in columns 11 and 9, 
respectively. 



cesses may freeze out at different radii, the simple 
picture of one energy sphere is certainly not en- 
tirely adequate. Therefore, it is quite surprising 
how well our simple description can account for 
the main features of the Monte Carlo results. 

3.3.2. Messer et al. (2001) Model 

As a second case we apply the same analysis 
to a model of Messer et al. (2001), provided to us 
by B. Messer. It is based on a full-scale Newto- 
nian collapse simulation of the Woosley & Weaver 
15 Mq progenitor model labelled sl5s7b, calcu- 
lated with the SN code described by Mezzacappa 
et al. (2001). We use a snapshot at 324 ms af- 
ter bounce when the shock wave is at a radius of 
about 120 km, i.e. the SN core still accretes mat- 
ter. We plot in Fig. 10 the radial profile of the 
density, medium temperature, temperatures T» 
and Tfj ux , particle flux, and luminosity (energy 
flux). Proceeding as before, we find that the en- 
ergy sphere should be at about 25 km (vertical 
dotted line) and that f (tes) ~ 40. The main con- 
clusions are the same as in the previous example. 

Our simple model for the and v T spectra ap- 
parently works so well because the predicted flux 
temperature is rather insensitive to the assumed 
location of the energy sphere. From the particle 
and energy fluxes in the bottom panel of Fig. 10 
one might have located the energy sphere, say, 
anywhere in the range tes = 23-29 km. One still 
would have predicted T nux within ±10% of the nu- 
merical value. Our approach of representing the 
energy sphere by a blackbody surface looks like a 
robust approximation. 

3.4. Is NN Bremsstrahlung Really Impor- 
tant? 

Within these specific models we can return to 
the question if nucleonic bremsstrahlung NN <-> 
NNvv is really important relative to the leptonic 
processes that were actually included in these sim- 
ulations. Based on the bremsstrahlung process 
alone we have calculated the thermalization depth 
Tthcrm as defined in Eq. (2), averaged over a ther- 
mal spectrum with the local medium temperature. 
For simplicity we have used the bremsstrahlung 
rate as in Eq. (C12) with B(x) = 1. We show 
T'thorm as the lower curve in the second panel of 
Figs. 9 and 10. 
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In the Janka & Hillebrandt (1989b) model 
Ttherm <C 1 at the location of the energy sphere so 
that bremsstrahlung would have been unimpor- 
tant even if it had been included. Bremsstrahlung 
depends on p 2 and thus drops quickly in low- 
density regions. This SN atmosphere has rather 
low densities so that it is not surprising that NN 
bremsstrahlung is unimportant. 

In the Messer et al. (2001) model the actual 
energy sphere is close to the bremsstrahlung ther- 
malization sphere, i.e. the radius where the lower 
curve in the second panel is 2/3. Therefore, in 
this model NN bremsstrahlung would have been 
roughly competitive with the leptonic processes 
that actually were included. Note that the den- 
sity at the energy sphere is about 10 times larger 
than in the Janka & Hillebrandt (1989b) case. 

In these models the atmospheres are still dom- 
inated by accretion. After a successful explosion 
they will settle to more compact dimensions, pre- 
sumably enhancing the role of NN bremsstrah- 
lung. We note that the model used by Hannes- 
tad & Raffelt (1998) represented a settled proto 
neutron star so that their finding of the impor- 
tance of bremsstrahlung is not inconsistent with 
the present discussion. Apparently it depends 
on the detailed atmospheric structure which pro- 
cesses predominate for the and v T thermaliza- 
tion. 

3.5. Sensitivity of Tfl ux to Thermalization 
Process 

We may explicitly test the sensitivity of the pre- 
dicted flux temperature Tfi ux to the strength of 
the assumed thermalization process by means of 
a power-law model of the stellar atmosphere as in 
Eq. (1), i.e. p oc r~ p and T oc r~ q . For the scat- 
tering and energy-exchange rates we assume 



A T oc p e cc r p e , 
A^ 1 oc (3p u T v t w cx/3r^ 



up — vq w 



(12) 



where for NN bremsstrahlung approximately u = 
2, v = 3 and w = —1. Further, [3 is a fudge 
factor which allows us to adjust the strength of 
the energy-exchange process. The thermalization 
rate then scales as 



(A^A^ 1 ) 



1/2 



withx = ^[(u+l)p+vq\ and y — ^(2+w). Finally, 
according to Eq. (10) we use 

Tflux cx T ES f- s * (14) 

with t = 0.128. 

It is now a matter of simple integrations to show 
that the transport optical depth and thermaliza- 
tion depth scale as r cx r - ^ -1 ) e 2 and T t horm oc 
pi/2 r -(x-i) e y rp^g i oca j thermal averages of the 
energy expressions are (e 2 ) cx T 2 cx r~ 2q and 
(e y ) cx T v oc r~ vq so that the thermally averaged 
optical depths are 

-(p+2 9 -l) 



Tithe 



CX 



oc 



-(x+yq-l) 



The energy sphere is located where r t he 
so that 



(15) 
= 2/3 



r ES OC ^(x+yq-1) ( 16 ) 

Therefore, the optical depth and temperature at 
the energy sphere scale as 

q 

T ES OC P 2(x+yq-l) ) 

p+2q-l 
TES OC (3 2(x+yq-l) _ 

Finally, with Eq. (14) we find 
T flux oc I3 Z 

with 



(17) 



(18) 



t{p + 2q-l)-q 



(19) 



2{x + yq-\) 
If NN bremsstrahlung is the thermalization pro- 
cess this translates into 

_ t(p + 2q-l)-q 
broms ~ 3p + 4q-2 ■ ( ' 

Since typically the density and temperature gra- 
dients are steep, dropping the —1 and —2 terms 
will introduce only a small error so that 

t(l + 2q/p) - q/p 



3 + 4q/p 



(21) 



oc f3 1/2 r~ x t v 



(13) 



Noting that typically we will have q/p — 0.25-0.35 
we find that Zbrcms is between —0.015 and —0.105. 
Taking this latter number as typical, a factor of 3 
uncertainty in the bremsstrahlung rate translates 
roughly into a 10% uncertainty of the predicted 
Tflux- Of course, this estimate depends on the tem- 
perature and density profile remaining unchanged 
by the variation of /?; a self-consistent treatment 
is not possible with the present approach. 



12 



3.6. Flavor Dependence of Spectral Tem- 
perature 

We may also address the question of the rel- 
ative spectral flux temperature between v e or v e 
and the other flavors for the case of a recoil-free 
scattering atmosphere. To this end we must spec- 
ify temperature and density profiles of the medium 
which we choose as power laws of the form Eq. (1) 
with T , r and po the energy-sphere values. For a 
plane-parallel geometry the medium temperature 
profile as a function of the transport optical depth 
r(e) for a given neutrino energy e is easily found 
to be 

T mod fl2Tj s T(e)\ u 



Ti 



ES 



e 2 T E s 



(22) 



with u = q/ (p — 1). If the density falls steeply so 
that p> 1 we have u ~ q/p. Realistic values for 
this power-law index are in the range, say, u = 
0.25-0.35. 

For a given energy e the v e or v e decouple 
at their respective energy spheres. The main 
energy-exchange process relevant for these species 
are charged-current (CC) processes involving elec- 
trons or positrons such as v e + p <-> n + e + . In 
addition there are neutral-current (NC) elastic col- 
lisions on nucleons. Since both CC and NC pro- 
cesses have the same e 2 behavior we may write 
^CC = £^nc w ith £ a factor which depends on the 
medium's composition and the degree of electron 
and nucleon degeneracy. Moreover, £ is different 
for v e and v e . In the absence of degeneracy effects 
we have £„ e « 4(1 — Y e ) and £ Pe « 4F e where 
Y e is the electron fraction per baryon. Note that 
the CC cross section is roughly 4 times the NC one 
and that only neutrons are available as CC targets 
for v e , and only protons for D e . We now apply 
the Shapiro- Teukolsky criterion Eq. (2) and find 
that the energy-dependent mfp for thcrmalization 

^ ^ therm = (ACC^NC) 172 = W + 0} ^ *NC ■ ^ 

tegrating over radius to obtain the optical depth 
then yields r th orm(e) = [£(1 + £)] 1/2r ( e ) where r(e) 
is the NC transport optical depth. The energy- 
dependent energy sphere is then located where 
7"thcrm( e ) = 2/3 so that the NC transport optical 
depth at this location is r(e) = (2/3)[£(l + £)r 1/2 - 
The above profile then informs us that 



Tj_ 

Tes 



(23) 



is the medium temperature at the energy depen- 
dent energy spheres for v e or v e . 

The main assumption for calculating the v e or 
v e spectra is that every energy group is emitted 
with a thermal flux corresponding to T e , i.e. we 
ignore chemical-potential effects. The escaping 
spectrum is then proportional to e 2 e -£ / Te . Note 
that every energy group is at the same transport 
optical depth because both NC and CC processes 
have the same energy dependence. (For every en- 
ergy group the energy sphere is at Ttherm(e) = 2/3 
so that this location also corresponds to a common 
value for the transport optical depth.) Therefore, 
the flux dilution caused by NC scatterings does 
not modulate the spectral shape. 

It is now a simple matter of integration to de- 
termine the v e or v e flux temperature and pinching 
parameter, 




irkr(f±|) 



r ( 4+2u\ 
1 V !+ 2 " J 



Vv e = 



( 4+2k\ p / 6+2tt 
^ l+2« J \ l+2« 



1 ^l+2u J 



(24) 



where T is the Gamma function. These spectra are 
always pinched because p < 1 for u > 0. Around 
the realistic value u = 0.3 we find the expansion 



0.19+Att 



Tflv, 



x (0.465- 0.65 Au), 
p„ e w 0.911 - 0.18 Au, 



(25) 



where Am = u — 0.3. 

As a first consequence we can estimate the rel- 
ative spectral temperature between v e and v e , 



Tfiu 



4F e (4T e + 1) 



4(i-y e )[4(i-y e ) + i] 



2(l + 2u) 



(26) 

For u = 0.3 and to lowest order in Y e this is 
0.69 (y/0.1) 094 , in the ball park of what is found 
in numerical simulations. 
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With the power-law representation of the 
flux temperature of Eq. (10) we find 



Tfi V 



(l.39 + 1.9 Auj 

_|_ ^jO.094+0.2 Au 
/f ES \ 0-06+0.4 Au 

\lo) 



(27) 



For v e we have £(1 + £) = 4F e (l + 4Y e ). Thus for 
small Y P we have to lowest order 



32+ 1.8 At* 



x (r e /o.i)°- 094+0 - 2Atl 

0.06+0.4 An 

X 



/ Tes \ 
V 10 J 



(28) 



This ratio is almost independent of the optical 
depth of the energy sphere; changing tes by 
a factor of 2 modifies the temperature ratio by 
around 5%. A change Am = 0.05 implies a 7% 
modification of the temperature ratio. Finally, 
changing Y e by a factor of 2 causes something like 
a 7% modification. 

In deriving this result we have made several 
simplifying assumptions so that the final answer 
should not be overinterpreted. Still, it is grat- 
ifying that our simple model, without any fine- 
tuning, reproduces relative neutrino flux tempera- 
tures which fully concur with what is found in nu- 
merical simulations. This finding, again, encour- 
ages us to put some faith in our explanation of the 
SN neutrino spectra. 

4. Nucleon Recoils 
4.1. First Example 

We now turn to the main goal of our investiga- 
tion, the effect of including nucleon recoils in the 
scattering atmosphere. To illustrate the type and 
magnitude of effect to be expected we begin with 
the same plane-parallel power-law model that was 
shown in Fig. 4, except that now we allow for recoil 
energy transfers in vN collisions according to the 
formalism described in Appendix B. The resulting 
radial runs of parameters are shown in Fig. 1 1 

Nucleon recoils cause the neutrino temperature 
T* to follow more closely the medium; the no- 
recoil case is shown as a long-dashed line. The 
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Fig. 11. — Example of Fig. 4 with nucleon re- 
coils included. The thin dashed lines refer to 
the case with a blackbody boundary condition 
Tes = 16 MeV at the energy sphere. The long- 
dashed line is T* for the no-recoil case of Fig. 4. 



Tflux profile is shifted downward so that the neu- 
trinos emerge with a lower Tfl ux at the surface. 
The luminosity profile has a "bump" near the en- 
ergy sphere, indicating that the neutrinos transfer 
energy to the medium within the scattering atmo- 
sphere. 

The energy sphere shown here is the same as 
in Fig. 11, i.e. strictly speaking we mean the ther- 
malization depth of NN bremsstrahlung when we 
say "energy sphere." Of course, in the present sit- 
uation the concept of an energy sphere no longer 
applies in the sense that nucleon recoils allow for 
energy transfers throughout the scattering atmo- 
sphere, which in turn no longer is strictly a scat- 
tering atmosphere. 

Next, we switch off NN bremsstrahlung en- 
tirely and use a blackbody boundary condition 
with T E s = 16 MeV at r E s- The resulting pro- 
files are shown as short-dashed thin lines. We 
practically obtain the same T* and T flux profiles 
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whether we allow the neutrinos to thcrmalize by 
NN bremsstrahlung or if we enforce a blackbody 
surface at tes , confirming our previous arguments. 

This approach has one more advantage. In 
general, ATfl ux depends on both the density and 
temperature profiles of the scattering atmosphere. 
However, if vN scattering is the only remaining 
process, then the only natural measure of distance 
is the transport optical depth r. Put another way, 
the medium density disappears from the Boltz- 
mann collision equation if we use a quantity pro- 
portional to JdrnB(r) as a radial coordinate. For 
us it is most practical to use the transport optical 
depth, averaged over a thermal neutrino spectrum 
at Tes- Note that we use the same thermal aver- 
age everywhere so that the T mec i profile does not 
enter the transformation from r to f . Then the 
only relevant profile is T me( j(f). 

Therefore, radial profiles such as Fig. 11 are ac- 
tually plotted in the "wrong" coordinate for issues 
related to neutrino transport. For illustration we 
show the same profiles in Fig. 12 with f as the ra- 
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Fig. 12. — Same as Fig. 11, except with the optical 
depth as the radial coordinate. 



dial coordinate, scaled to tes- We recognize that 
all temperature profiles vary slowly in the neigh- 
borhood of the energy sphere. 

The "bump" in the luminosity profile turns 
out to extend throughout the scattering atmo- 
sphere. For the blackbody boundary condition, 
the luminosity decreases nearly linearly toward the 
surface so that the neutrinos transfer energy to 
the medium throughout the scattering atmosphere 
with approximately the same rate everywhere. It 
is interesting that the luminosity profiles are very 
different between the 7V./V-tliermalization case and 
that with a blackbody boundary condition. 

4.2. Constant Temperature 

Turning to a systematic study of the impact 
of nucleon recoils we consider a scattering atmo- 
sphere with a fixed temperature T me( j, whereas the 
blackbody surface at the bottom is at the temper- 
ature Tes > T me d. What is the spectrum of the 
emerging neutrinos? 

Actually not even the sign of ATfl ux caused by 
nucleon recoils is a priori obvious. Since Tg ux < 
Tes it is conceivable that energy exchange will in- 
crease Tfl ux in a situation when T me( j > Tg ux . (The 
superscript refers to the no-recoil case.) Con- 
versely, one expects the energies of the emerging 
neutrinos to be lowered when T mod < T^ ux . 

The real answer is more complicated. In Fig. 13 
we show the T* and Tfl ux profiles for a scattering 
atmosphere with f E s = 100 and T ES = 20 MeV. 
The dashed lines are for the no-recoil case. The 
three panels are for the different medium temper- 
atures 20, 12, and 4 MeV as indicated by dotted 
lines. 

In the first panel where T mc d = Tes the trapped 
neutrinos approach more closely T moc i, while Tfl ux 
at the surface actually increases! In the second 
panel we use T^ ux < T mod < T E s- The T* profile is 
close to T moc i while Tfl ux shoots up near the energy 
sphere. However, in spite of T^ ux < T mo d, at the 
surface Tfl ux reaches a value below the no-recoil 
case. The overall behavior is similar in the third 
panel where T mod < T,? ux < T ES . 

In Fig. 14 we show the T» and Tfl ux profiles 
for T mcd varying from T ES = 20 MeV to 2 MeV 
in steps of 2 MeV. While the T» curves certainly 
look vaguely as expected, the qualitative behavior 
of the Tfl ux profiles would not have been intuitively 
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Fig. 13. — Neutrino temperature profiles T* and 
T flux for T ES = 20 MeV and T mod = 20, 12 and 
4 MeV (top to bottom). Dashed lines for the no- 
recoil case, dotted lines for T mc d. 



obvious. Interestingly, the T nux curves all intersect 
essentially in one point. 

Finally we show in Fig. 15 for several values 
of tes how the parameters of the emerging neu- 
trino flux vary with T mc d. In the first panel we 
show Tft ux , in the second panel the relative change 
ATflux/Tfl^ where the superscript again refers to 
the no-recoil case. These latter curves intersect es- 
sentially in one point, suggesting that the medium 
temperature for which the neutrino flux spectrum 
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Fig. 14. — Neutrino temperature profiles T* and 
T flu x for T ES = 20 MeV. From top to bottom the 
curves are for T mo d = 20 to 2 MeV in steps of 
2 MeV. Dashed lines arc for the no-recoil case. 



remains unchanged is nearly independent of the 
optical depth of the scattering atmosphere. 

In the third panel we show the pinching param- 
eter. The effect of nucleon recoils is not only to 
shift Ta U x, but also to modify the relative shape of 
the spectrum. If the difference between Tes and 
T me d is large enough the spectra become pinched. 

In summary, we have found that nucleon recoils 
can have the effect of shifting Tfl ux of the emerg- 
ing neutrinos both up or down, depending on how 
much lower T nU x is relative to T mo d- Likewise, the 
spectral shape can become more or less pinched. 

In a SN core the temperature varies slowly as a 
function of f, i.e. for most of the scattering atmo- 
sphere T mo d is close to Tes (Fig. 12). Therefore, it 
is by no means obvious if nucleon recoils will shift 
the energies of the emerging neutrinos up or down, 
even though naively one would have expected that 
an energy-transfer channel between neutrinos and 
nucleons must always shift them downward. A 
simple inspection of the temperature profiles in 
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Fig. 15. — Variation of neutrino flux parameters at 
the surface with medium temperature for several 
values of tes- 



the no-recoil case does not allow one to predict 
the direction of the energy shift. This unintuitive 
behavior is a consequence of the huge difference 
between the neutrino temperature and their flux 
temperature, which in turn is a consequence of the 
energy dependence of the vN scattering cross sec- 
tion. If this cross section were energy independent, 
these complications would not arise. 



4.3. Power-Law Profiles 

As a next step we consider temperature profiles 
which resemble a realistic SN core, i.e. power-laws 
of the form 



- mod 



(29) 



If the density varies with radius as p cx r~ v and 
T cx r~ 9 , then u = q/(p — 1). The density is 
a steeply decreasing function of r, i.e. p » 1 so 
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that u w q/p. A realistic range is q/p = \ - 
i, i.e. the relative temperature profile does not 
depend sensitively on p. In the following we will 
consider the range < u < 1 where the extreme 
case u = corresponds to a constant T moc i = Tes 
while the other extreme u = 1 corresponds to T me <i 
decreasing linearly from the energy sphere to the 
surface. 

For tes = 30 we show in Fig. 16 profiles for 
Tmed, T, and Tfl ux , from top to bottom for u = 0, 
0.3 and 1. We also show as dotted horizontal lines 
the no-recoil value Tg ux . In Fig. 17 we show the 
Tflu X profiles for u = to 1 in steps of 0.1. For 
u > 0.1 the neutrinos emerging from the surface 
have a lower temperature than in the no-recoil 
case. For a realistic power-law index u m 0.3, Tfl ux 
decreases nearly linearly from the energy sphere 
to the surface. It deserves mention that relative 
to the no-recoil case Tfl ux near the energy sphere 
has increased, at the surface it has decreased. 
For power-law indices which arc not too small, 
the T nux profiles roughly intersect in one point, 
roughly halfway between the energy sphere and 
the surface. 

In Fig. 18 we show the relative shift AT Rux /T^ ux 
of the surface flux temperature caused by nucleon 
recoils as a function of u for tes = 10, 30 and 100. 
The shift is upward for u < 0.1, the exact value 
depending on tes- For realistic values around 
u = 0.3 the shift is between about —5% and —20%, 
depending on tes- We also show the pinching pa- 
rameter which is larger than 1 for small u, indicat- 
ing anti-pinched spectra. For large u the spectra 
are pinched. In the realistic range around u = 0.3 
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Fig. 17. — Tfiu X profiles for u = to 1 in steps of 
0.1. 
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Fig. 18. — Spectral properties at the surface for 
power-law models. 

the pinching parameter is slightly below 1 so that 
the spectra are pinched, but rather close to a ther- 
mal shape. 

4.4. Variation with Nucleon Mass 

Next we address the question of how sensitive 
ATflux depends on the efficiency of recoil energy 
transfer. The only dimensionless parameter gov- 
erning recoil is T/m with m the nucleon mass. 
Therefore, we may either change the absolute scale 
of the medium's temperature profile, or we may 
artificially change the nucleon mass. 

We have opted for varying m and show in 
Fig. 19 profiles for T* and Tfl ux for a constant 
T mcd = 5 MeV and for f ES = 100. Short- 
dashed lines show the no-recoil case, correspond- 
ing to m = oo, while long-dashed lines refer to 
the standard recoil case with m = mo where 
too = 938 MeV is the physical nucleon mass. The 
solid curves are from top to bottom for nucleon 
masses to/to = 10 4 , 10 3 , 10 2 , 10 and 10" 1 . 

Surprisingly, even a huge to allows for a large 
modification of the local neutrino temperature. 
Neutrinos scatter frequently so that even small 
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Fig. 19. — Neutrino temperature profiles T* and 
T flux for T ES = 20 MeV and T mcd = 5 MeV 
(thin dotted line). From top to bottom the curves 
are for nucleon masses to = oo (no-recoil, short 
dashes), m/m = 10 4 , 10 3 , 10 2 , 10, 1 (standard 
case, long dashes), and 10 _1 . 
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Fig. 20. — Spectral properties of escaping neu- 
trinos for constant T mcc j as indicated, and with 
T ES = 20 MeV. 



energy transfers compound to a large effect. In 
Fig. 2 we showed the spread in final-state energies 
for a typical neutrino interacting with thermal nu- 
cleons. It is worth noting that the width of this 
curve scales with (T/to) 1 / 2 so that the analogous 
curve for to = 10 4 too would have 1% of the origi- 
nal width. 

This is further illustrated in Fig. 20 where 
we show the m-dependent Tfl ux of the escaping 
neutrinos for f E s = 100, T ES = 20 MeV and 
Tmcd = 5, 10 and 20 MeV. For m -> oo (no- 
recoil) the asymptotic value is 9.30 MeV. When 
Tmed = 20 MeV recoils shift Tfl ux upward, for 
Tmcd = 5 MeV downward, both effects being 
monotonic with to. For T mo d = 10 MeV the shift 
is downward, but Tfl ux is not a monotonic function 
of m. In the lower panel we show the variation of 
the pinching parameter with to. 

In Fig. 21 we perform the same analysis for a 



power-law medium with u = 0.3, Tes = 20 MeV, 
and tes = 100. Again, Tfl ux of the escaping neu- 
trinos is not a monotonic function of to; the shift 
ATg ux relative to the no-recoil case is most ex- 
treme for to around the physical value. Moreover, 
the pinching parameter shown in the lower panel 
is quite close to 1 for to near its physical value 
so that the spectrum takes on a nearly thermal 
shape. We believe it is coincidental that a nucleon 
mass near its physical value maximizes the energy- 
transfer effect for the conditions of interest. 

In summary, the energy shift of the escaping 
neutrinos and their spectral shape are both ex- 
tremely insensitive to the nucleon mass, i.e. to the 
exact width of the medium's dynamical structure 
function for energy exchange. Changing the nu- 
cleon mass by a factor of a few up or down has 
no significant effect whatsoever on the predicted 
neutrino flux spectrum. 
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Fig. 21. — Spectral properties of escaping neutri- 
nos for a medium with a power-law T mo d profile 
with u = 0.3, T ES = 20 MeV, and f ES = 100. 



4.5. Is Inelastic uN scattering important? 

Many-body effects modify the nuclear medium's 
dynamical structure function. However, even 
when one ignores nucleon-nuclcon correlations, 
neutrinos can transfer energy by "inelastic scat- 
tering" vNN — > NNv 1 the crossed version of 
bremsstrahlung NN —> NNvv. This channel al- 
lows for energy transfer even if one ignores recoil. 
For this effect one may use the same dynamical 
structure function as for bremsstrahlung — see Ap- 
pendices C and D for technical matters. Brems- 
strahlung varies as p 2 so that it depends on the 
density if recoil alone or inelastic scattering alone 
is more important. 

In view of the insensitivity of Tfl ux to the exact 
width of the structure function it is clear that in- 
elastic scattering can not have a large impact in 
addition to recoil. We have actually performed a 
variety of runs where we include the combined ef- 
fect of recoil and inelastic scattering by means of 
the combined structure function described in Ap- 
pendix D. For example, in a power-law model like 



the one shown in Figs. 4 and 11, including only the 
inelastic vN channel gives a significant but weaker 
shift of Tflux than recoil alone. Combining both ef- 
fects, Tfl ux is virtually indistinguishable from the 
recoil-only case. 

Again, we conclude that the detailed implemen- 
tation of energy transfer does not matter as long 
as its efficiency is crudely comparable to that of 
recoil. 

5. Discussion and Summary 

We have studied how nucleon recoils in a SN 
core's scattering atmosphere modify the escaping 

and y T spectra. We justified and used a model 
where the energy sphere at the bottom of the scat- 
tering atmosphere is represented by a blackbody 
boundary condition for the neutrino distribution 
function (temperature Tes)- The only remaining 
interaction process is vN scattering. If these col- 
lisions are taken to be iso-energetic, the proper- 
ties of the emerging neutrino flux are determined 
by exactly two parameters, Tes and the transport 
optical depth at the energy sphere tes, averaged 
over a thermal spectrum. For spherically symmet- 
ric rather than plane-parallel geometry, the analo- 
gous role of f E s is played by a quantity where the 
column density is calculated by including a factor 
(rEsA) 2 under the integral (energy-sphere radius 
r E s)- However, for the steep density gradients typ- 
ical for SN cores the difference between a spher- 
ical and plane-parallel treatment is small. Real- 
istic SN models have tes ~ 10-30. In this range 
the spectral flux temperature Tfl ux of the escaping 
neutrinos is 62-54% of Tes- Equation (9) is an 
analytic approximation formula for Ta ux /Tes as a 
function of tes- In summary, the spectral temper- 
ature of the and v T flux emitted from a SN core 
is roughly 60% of the medium temperature at the 

and v r thermalization depth. This is confirmed 
by SN simulations from the literature. 

The effect of nucleon recoils on Tg ux is in sev- 
eral ways counter-intuitive. We always assume 
monotonically decreasing temperature profiles so 
that throughout the scattering atmosphere T moc i < 
Tes- Still, the average emerging neutrino ener- 
gies can be both larger or smaller compared with 
the no-recoil case. This behavior is explained 
by the large difference between the local tem- 
perature T* of trapped neutrinos and their spec- 
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tral flux temperature Tfl ux , a difference which is 
caused by the e 2 dependence of the vN scatter- 
ing cross section. If we assume a power-law model 
T-med = Tes(t/tes) u > the spectral shift is upward 
for u < 0.1 (Fig. 18). Realistic power-law indices 
are around u = 0.3 where the shift is downward 
by anywhere between —5% and —20%, depending 
on u and tes- 

The spectral shift is extraordinarily insensitive 
to the assumed value of the nucleon mass m, i.e. 
to the details of the energy-transfer mechanism. 
Therefore, for the spectra formation problem we 
do not need to know the details of the nuclear 
medium's dynamical structure function. As a con- 
sequence, the relative shift of Tfl ux is also insensi- 
tive to the absolute scale of T me( j, it is only sensi- 
tive to its radial profile, in our case the power-law 
index u. For power-law profiles, Tfl ux /T E s depends 
on tes and u, but only extremely weakly on T E s 
or m. 

The magnitude of the spectral shift found in 
our investigation is too large to ignore in full-scale 
numerical simulations, but too small to be wildly 
critical. If the relative flux spectra, say, between 
v e and differ by only 10-30% as in some sim- 
ulations, then a reduction of the v^ energies by 
10-20% would be quite important. If, on the other 
hand, the spectra are much more different as found 
in other simulations, then the recoil reduction may 
not make much of a practical difference. The im- 
portance of recoil effects may well depend on the 
exact phase within the SN evolution. There may 
be big differences between the post-bounce pre- 
cxplosion accretion phase, and the later Kelvin- 
Helmholtz cooling phase. 

We also warn that our discussion, being based 
on a static background model without self- 
consistent adjustment, may be more relevant as 
an upper limit to the possible magnitude of the ef- 
fect than a realistic estimate. The energy transfer 
permitted by nucleon recoils allows the neutri- 
nos to heat the medium, modifying T mc d such as 
to reduce the rate of energy transfer. Of course, 
such a modified temperature profile will affect 
the electron-flavored neutrinos, perhaps increas- 
ing their energies. We expect that nucleon recoils 
cause the v e and v^ spectra to be more similar 
than they would otherwise be, but it would not 
be correct to estimate the multi-flavor spectra by 
taking the spectra from a simulation and reducing 



the v^i energies by our recoil factors. 

Another problem with the interpretation of our 
result is that we have always neglected ve scatter- 
ing. This process provides a large amount of en- 
ergy exchange in a collision and thus should freeze 
out at its thermalization sphere. Nucleon recoils, 
on the other hand, provide for a small amount of 
energy transfer in each collision, a process which 
does not freeze out, except when neutrinos stream 
freely beyond the transport sphere. Therefore, one 
would expect the role of vN and ve scattering to 
be rather different. We expect our treatment to 
apply beyond the ve freeze-out sphere. However, 
given the extreme insensitivity of the spectral shift 
to the exact mode of energy transfer, we feel not 
entirely sure about what happens when ve scatter- 
ing and vN recoils are both present, even if the ve 
freeze-out sphere is in the neighborhood of the NN 
bremsstrahlung sphere. The effect of nucleon re- 
coils has turned out to be rather counter-intuitive 
in several ways, cautioning us not to rush into pre- 
dictions about the behavior of this coupled system 
without having investigated it. Again, it appears 
safe to take our estimates as upper limits to the 
differential change expected when nucleon recoils 
are included in addition to the ve process. 

In a nutshell, the main conclusions of our in- 
vestigation are that the v^ and v T spectra are well 
accounted for by the picture of a scattering at- 
mosphere with a blackbody boundary condition 
at the bottom, that nucleon recoils lower the flux 
temperature for typical conditions by as much as 
—20% if all else is kept equal, and that this effect is 
astonishingly insensitive to the detailed treatment 
of the energy-transfer mechanism. 
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A. Neutrino Transport in the Diffusion 
Limit 

A.l. Plane Parallel Geometry 

In order to calculate the flux suppression, or 
"flux dilution," by a scattering atmosphere in the 
diffusion limit we assume that the neutrinos inter- 
act with "infinitely heavy" (recoil free) scattering 
centers. The cross section is taken to be 



da 
dcos 9 



Co 



1 + &cos 6 



(Al) 



where — 1<6<+I,0is the scattering angle, and 
cto the total scattering cross section. The neutrino 
energy e is assumed to be fixed so that the energy 
dependence of <j a can be ignored. 

We ask for the neutrino flux penetrating a cer- 
tain medium layer. To this end we need to solve 
the Boltzmann Collision Equation for the distri- 
bution function / of the neutrinos, 



8f 



V/ = C[/], 



(A2) 



where v is the neutrino velocity and C[f] the col- 
lision term. When the neutrinos do not exchange 
energy with the medium this equation applies to 
each neutrino energy e separately so that indeed 
we may assume e to be fixed. 

As a first case we assume that the medium is 
layered plane-parallel transverse to some direction 
r. In a stationary state and taking massless neu- 
trinos (|v| = 1) the collision equation is 

pd r f(r,p) = -a n B (r)f(r,p) 

r+l 

+ n B {r) j dp' {a)^f{r,p!), 

(A3) 

where n B (r) is the number density of scattering 
centers ("baryons"). The distribution function de- 
pends only on the one spatial coordinate r — |r| 
and on p, the cosine of the propagation angle rel- 
ative to r. The azimuthally averaged differential 



scattering cross section (a) ^ represents the prob- 
ability for a neutrino to scatter from the direction 
p! to p. If the cross section is given by the "dipolc 
formula" Eq. (Al) we find (a)^ = a (l + bpp')/2 
so that the r.h.s. of the collision equation becomes 



a Q n B (r) 



r +i 1 



bpp! 



(A4) 



We may further use the spatial coordinate 

p=a Q [ dr'n B (r'), (A5) 

J r a 

corresponding to the optical depth at location r as 
measured from the radiating surface at the bottom 
of the medium layer (r = ro). The collision equa- 
tion then simplifies to 

fJ-d p f(p^) = -f(p,p) 

+1 , 1 + bpp' , 

dp ^ /(P,M)- 

-l 1 

(A6) 

To solve it, we must supplement it with the bound- 
ary conditions f(0,p) = <f>o for < p < 1 and 
f(r,p) — for — 1 < p < at the medium's sur- 
face which is characterized by r = p(oo), the total 
optical depth of the scattering atmosphere. 

An exact analytic solution of Eq. (A6) is found 
by an ansatz corresponding to the usual moment 
expansion 



f(p,p) = $ 



f 



P 

T 



a(l + p) , (A7) 



where 



a- 1 =r (1-6/3). (A8) 

If the medium is optically thick, t>1, the bound- 
ary conditions are approximately fulfilled. 

The first moment of the distribution func- 
tion, which is proportional to the particle flux, is 
27r JJ^ 1 dp p f(p, p) and is found to be (47r/3) a $o, 
independently of p, i.e. the flux is conserved. If 
the medium were absent, the solution would be 
f(p, p) = $o©(m) with the flux 7r$ - Hence in the 
diffusion limit we find 



s(t) 



(3-6)r 



(A9) 



for the neutrino flux suppression relative to the 
free-streaming case. 
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We observe that the r.h.s. of Eq. (A8) corre- 
sponds to the usual transport cross section, 



Therefore, the spherically symmetric solution is 
obtained by the substitution 



r+l 

C7 T = / d cos 9 (I — cos 9) 



da 
dcos9 



(A10) 



With Eq. (Al) this is a T = cr (l - b/3). There- 
fore, the relevant variable is the transport cross 
section, not the scattering cross section. In the 
main text we usually mean ctt when we speak of 
the cross section, and we mean t(1 — b/3), the 
"transport optical depth," when we speak of the 
optical depth. 

A. 2. Spherical Symmetry 

Next, we turn to the case of spherical symmetry 
instead of the previous plane-parallel geometry. If 
one writes the neutrino distribution function in 
terms of the variables r and /i, where [i is the 
cosine of the propagation angle relative to the ra- 
dial direction, the l.h.s. of the stationary collision 
equation (A3) becomes 



1 



f(r,fi). 



(All) 



With the moment expansion f(r, /x) = fo{r) 
H fi (r) the collision equation turns into 



Hd r f + n 2 d r f 1 + 



1 



/l = -GTTlBllh- (A12) 



<J T n B (r) -» a T n B (r)(r Q /r) 2 



(A15) 



Put another way, in the diffusion limit the scat- 
tering atmosphere dilutes the neutrino flux by the 
factor 

«(r.) = ^-, (A16) 
3r» 

where 

r* = a T ^ drn B (r)(^y. (A17) 

Of course, even in the spherical case the proper 
optical depth is defined without the (r /r) 2 factor 
under the integral. However, for the purpose of 
neutrino flux transmission the parameter t* plays 
the same role as the proper optical depth r in the 
plane-parallel case. 

We may compare the two cases for an example 
where the density drops as a power law, n B {r) cx 
r~ p , implying 

t* p-l 



p+1 



(A18) 



If p » 1 this ratio is close to unity. Put another 
way, if the density drops quickly for r > r the 
integral is dominated by r-values near r and we 
recover the transport optical depth of the plane- 
parallel case. 



The neutrino flux is (47r/3) /i(r) so that flux con- 
servation through a given spherical shell implies 
/i(r) cx r~ 2 or 



fJ*d r fo(r) 



1-3 M 2 



fi(r). 



(A13) 

In the diffusion limit, the second term in square 
brackets is much larger than the first for all r 
which significantly exceed the scattering mean free 
path. Neglecting the first term we find 



d r fo(r) = -<T T n s (r)/i(r). 



(A14) 



The only remaining difference to the plane-parallel 
case is that flux conservation now dictates that 
fi(r) = <&oa(ro/r) 2 with ro the bottom of the 
scattering atmosphere. In the plane-parallel case 
we had /i(r) = $o a. 



B. Neutrino-Nucleon Scattering 

B.l. Scattering Cross Section 

We need to derive an expression for the dy- 
namical structure function S(u>, k) which plays the 
role of a scattering kernel for neutrino processes. 
In terms of the structure function the differential 
scattering cross section of a neutrino with initial 
energy ei is 



da 



de 2 d cos 9 



Ci(3-cos0) 2 2 S(u,k) 
- (_r F e 2 



2tt 



2tt 



(Bl) 



where e 2 is the final-state neutrino energy, lo = 
ei — e 2 the energy transfer, k = |ki — k 2 | the mod- 
ulus of the momentum transfer, and 9 the scatter- 
ing angle. S(u, k) is a function of k, and not of k, 
because the medium is assumed to be isotropic. 
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In this paper we include only one species of 
nucleons for which the axial weak coupling con- 
stant is assumed to be \Ca\ — 1.26/2. In general, 
there is also a vector-current interaction; the total 
scattering cross section on non-relativistic nucle- 
ons is proportional to Cy + iC\. For protons, 
Cy very nearly vanishes while even for neutrons 
with Cy = —1/2 the vector-current contribution 
to the cross section is a rather small correction. 
Ignoring the vector current simplifies our discus- 
sion as we need only one structure function. In 
general, there is a different structure function for 
the vector-current, the axial vector current, and 
the mixed term. 

When the nucleons are "infinitely heavy" and 
thus unable to recoil, the structure function is sim- 
ply S(oj,k) = 2ttS(uj), i.e. the collisions are iso- 
energetic (no energy exchange), and the nucleon 
can absorb any amount of momentum. In this case 
the total axial-current scattering cross section and 
the transport cross section for a neutrino of energy 
e are 



iC 



(7t = 



— (_rp£ , 



7T 

3tt f 



(B2) 



This corresponds to an inverse mean free path of 



12.08 km 1 pue\ , 



Ao 1 

At 1 = 13.42 km" 1 p 14 e 2 w , (B3) 
where pu = p/10 14 g cm~ 3 and e w = e/10 MeV. 

B.2. General Properties of the Structure 
Function 

The dynamical structure function can be ex- 
pressed as a current-current correlator, in our case 
of the nucleon axial vector current (see, for exam- 
ple, Janka et al. 1996). This representation al- 
lows one to derive a number of general properties. 
First, there is detailed balance 



S(-u,k) 



-ui/T 



S(w,k). 



(B4) 



Second, if there are no spin-spin interactions, we 
have the normalization 



^-S(w,k) = l, (B5) 



and the /-sum rule 



S 

J — ( 



lu S(lu, k) = - — . 
2tt x 1 ' 2m 



(B6) 



Recent calculations of the dynamical structure 
functions in the presence of nucleon-nucleon in- 
teractions, i.e. in the presence of spin-spin cor- 
relations, have been performed by Burrows and 
Sawyer (1998) as well as Rcddy, Prakash and Lat- 
timer (1998) and Reddy et al. (1999). In our 
present study we ignore all nucleon-nucleon inter- 
action effects except for the calculation of nucleon- 
nucleon bremsstrahlung (sec below). 

It is sometimes useful to construct a symmetric 
version of the structure function, 



1 + e-"/ T 



S(u,k), (B7) 



and conversely 

S(cu,k) 



1 + e-"/ T 



S(w,k). 



(B8) 



By construction S(oj, k) is a symmetric function of 
uj; it obeys the same normalization as S(ui, k). 

B.3. Nucleon Recoil 

For a finite nucleon mass and ignoring nucleon 
degeneracy effects, the structure function takes on 
the well-known form 



•Srecoil^, fc) — 



UJ k T 



exp 



4Tw fc 



, (B9) 



where uJk = k 2 /2m. One easily checks that 
this "displaced Gaussian" fulfills detailed balanc- 
ing Eq. (B4) as well as the sum rules Eqs. (B5) 
and (B6). 

The distribution da/de.2 OI final-state neutrino 
energies is obtained by integrating Eq. (Bl) over 
cos# with Eq. (B9) for S(w, k) and observing that 
lu = e 1 — e 2 and k = (e\ + e 2 — 2eie 2 cos ff) 1 ! 2 . 
As an example we show in Fig. 2 this distribution 
for T = 10 MeV and an initial-state neutrino en- 
ergy ei = 30 MeV. Evidently the distribution is 
rather broad. For other values of t\ it looks simi- 
lar when the horizontal axis is scaled accordingly, 
i.e. the fractional width relative to ei does not de- 
pend much on e\. 
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C. Bremsstrahlung 

C.l. Nucleon Spin Relaxation Rate 

In order to calculate the bremsstrahlung emis- 
sion of neutrino pairs NN — > NNvv and re- 
lated processes we need to include nucleon-nucleon 
interactions. The neutrino energy-loss rate of 
a single-species non-relativistic non-degenerate 
thermal nucleon gas can be expressed as (see, e.g., 
Raffelt 1996) 



d 3 ki 



n B 



d 3 k 2 



(ei + e 2 ) 



2ei(27r) 3 2e 2 (27r) 3 
x 8eie 2 (3 - cosfl) S (-ei-e 2 , |ki+k 2 |) , 

(CI) 

where ei i2 are the energies of the emitted neutri- 
nos, ki, 2 their momenta, and S(u>, k) is the same 
dynamical structure function that describes axial- 
current neutrino-nucleon scattering. In the non- 
relativistic limit the vector current does not con- 
tribute to bremsstrahlung. 

The phase-space integration in Eq. (CI) cov- 
ers energy-momentum transfers (u>, k) = (— ei— e 2 , 
|ki+k 2 |) which are time-like (w 2 > k 2 ) whereas 
for neutrino scattering they are space-like (ui 2 < 
k 2 ). Therefore, even though both scattering and 
bremsstrahlung are characterized by the same 
S(uj,k), it is different regions in the (w,fc) plane 
that contribute to these processes. The recoil 
structure function Eq. (B9) has only power for 
space-like (u), k) and thus cannot account for 
bremsstrahlung, in keeping with the obvious in- 
sight that free nucleons cannot radiate. Nucleon- 
nucleon interactions were included in various cal- 
culations of S{oj, fc), most recently by Burrows and 
Sawyer (1998) as well as Reddy, Prakash and Lat- 
timer (1998) and Reddy et al. (1999). However, 
these and previous works limited their calcula- 
tions to the random-phase approximation where 
nucleon spin-spin correlations are included, but 
not nucleon spin fluctuations. Bremsstrahlung 
requires nucleon spins to fluctuate. Graphically 
speaking, a nucleon spin needs to be kicked to be 
able to emit radiation which couples to the spin 
such as neutrino pairs. The S(uj 7 k) calculated in 
these works have no power for time-like (w, k) and 
thus do not account for bremsstrahlung. 



One finds essentially two different approaches 
in the literature to calculating the bremsstrahlung 
rate. One may model the nucleon-nucleon inter- 
action potential, typically by a one-pion exchange 
potential, and then evaluate the relevant Feynman 
amplitudes for NN —* NNvv. Recently, Hanhart, 
Phillips and Reddy (2001) have followed a different 
approach where they calculate the bremsstrahlung 
matrix element from nucleon scattering data in the 
soft limit (uj — > 0) and then extend the results "by 
hand" to non-vanishing energy transfers. 

If one follows this latter approach one needs to 
guess a functional form for S(u,k). To this end 
one may use the "long wavelength approximation" 
where the momentum transfer k to the nucleons is 
assumed to be negligible so that one approximates 
the structure function by S(u) = lim fe ^ S(oj, k). 
We further recall that the structure function can 
be expressed as a nucleon spin-spin correlator. 
The nucleon spin autocorrelation function has a 
nontrivial time dependence only if there are spin- 
nonconserving forces between the nucleons such as 
the nuclear tensor force. It is plausible (but not 
necessary) that a nucleon spin, being kicked by 
other nucleons, loses memory of its original ori- 
entation as e _r * where T is the spin relaxation 
rate. An exponential decline corresponds to the 
assumption that the spin relaxation is a Marko- 
vian process, i.e. a chain of uncorrelated random 
kicks. The Fourier transform of an exponentially 
declining autocorrelation function is 



2T 



+ r 2 ' 



(C2) 



leading to a Lorentzian ansatz for the symmetric 
structure function. 2 

Calculating the bremsstrahlung rate in the one- 
pion exchange model in Born approximation with 
uncorrelated, non-degenerate, single-species nucle- 
ons, Raffelt & Seckel (1995) find in the soft limit 
(u> — > 0) and ignoring the pion mass 



T = 2y / 7ra 2 hbT 



V2 m -B/2 



(C3) 



2 In previous papers (Janka ct al. 1996 and references) we 
have used the notation r CT = 2T and called r CT the spin fluc- 
tuation rate. However, the relationship between V and an 
exponentially declining spin autocorrelation function (Raf- 
felt & Sigl 1999) implies that T, and not r CT , has the intu- 
itive interpretation of a spin relaxation or spin fluctuation 
rate. 
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with ub the baryon density, m the nucleoli mass, 
and = (/2m/TO7r) 2 /47r 15 with / ss 1 the 
pion-nucleon "fine-structure constant." 

In the one-pion approximation one can calcu- 
late the behavior of the structure function for non- 
vanishing energy transfers and one may also in- 
clude the pion mass. Therefore, Eq. (C2) can be 
supplemented with a dimcnsionless factor s(x, y) 
where x — ui/T and 



(C4) 



with Tio = T/10 McV. One finds (Raffelt & Seckel 
1995; Hannestad & Raffelt 1998) 



s(x,y) = 



-x/2 + e x/2 



/'OC 

J\x\ 



16 



dt £ _ t/2 3x 2 + 6ty + 5y 2 



2V^" 



x 2 + 2ty + y 2 
1 , ft + y + Vt 2 - x 2 



t + y 



log 



t + y- Vt 2 - 



(C5) 



This expression is even in x and also s(0, 0) = 1. 

In Fig. 22 we show s(x, y) as a function of 
x for several values of y. In the outer parts of 
the SN core, we have y = 2-4, values for which 
s(x, y) varies slowly as a function of x. The lower 
panel of Fig. 22 shows that in this y-range we have 
2y/T+ys{x,y) = 1 within about ±10%. There- 
fore, at the relatively crude level of approximation 
that we content ourselves with it is justified to use 
the Lorentzian Eq. (C2) as a structure function 
with the spin relaxation rate Eq. (C3), divided by 
2VTTy, i.e. 



^/Trains T 



m 2 y/ mT + m 2 
Within a few percent this is numerically 



7; 



= 1.25pi4 



2 + T, 



10 



(C6) 



(C7) 



where p w = p/10 14 g cm" 3 and T w = T/10 MeV. 

It is not obvious how well or poorly this re- 
sult represents the true spin relaxation rate for 




Fig. 22. — Dimensionless structure function ac- 
cording to Eq. (C5) for the indicated values of y. 



the physical conditions at hand. Hanhart, Phillips 
and Reddy (2001) have calculated the equivalent 
of T on the basis of nuclear scattering data. They 
claim that at nuclear density the true T is about 
a factor of 4 smaller than the one obtained from 
the OPE calculation. 3 In a one-species medium 
at nuclear density (3 x 10 14 g cm -3 ) the nucleon 
Fermi momentum is about 350 MeV, much larger 
than the pion mass, so that nucleon-nucleon col- 
lisions probe deeply into the core of the nucleon- 
nucleon potential; it is no surprise that the OPE 
model, which only mimics the peripheral part of 
the potential, would yield poor results. However, 
we are here interested in a non-degenerate medium 
at, say, T = 10 MeV where a typical nucleon mo- 
mentum is {ZmT) 1 / 2 w 170 MeV so that the OPE 



3 These authors compare with the OPE calculations of 
Friman and Maxwell (1979). However, Friman and 
Maxwell have missed a term in the squared matrix ele- 
ment of NN — > NNuP; the correct result was derived by 
Brinkmann and Turner (1988) as well as Raffelt and Seckel 
(1995). Therefore, the discrepancy between the OPE re- 
sult and that of Hanhart, Phillips and Reddy (2001) is even 
larger, perhaps by as much as a factor of 1.5. 
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method would be expected to be much more justi- 
fied. Unfortunately, Hanhart et al. do not provide 
an estimate for such conditions. Hannestad & Raf- 
felt (1998) have argued that the OPE result should 
not underestimate the true answer by much more 
than about 30% for such conditions. 

An additional complication arises in a mixed 
medium of protons and neutrons. Within the 
OPE model, the bremsstrahlung rate of such a 
mixed medium far exceeds that of a single-species 
medium at the same density because the np ma- 
trix element is much larger. In the present paper 
the medium is always modeled by a single-species 
composition, leading to an underestimate of the 
bremsstrahlung rate of a realistic SN core. It is 
not quantitatively known how well these two ef- 
fects compensate each other, but we would have 
to be very unlucky if our adopted value for the 
spin relaxation rate were off by more than a fac- 
tor of 2 in either direction. 

C.2. Energy Loss Rate 

We may now calculate the energy-loss rate of 
a medium due to nucleon-nucleon bremsstrahlung 
emission of neutrino pairs. Equation (CI) can be 
simplified to read 

Q v v = -^n B du^Si-uj), (C8) 



where the Lorentzian structure function is 
2r 2 



S(w) 



(C9) 



T 2 +oj 2 i + e --v T ' 

Since S(ui) is multiplied with a high power of u 
under the integral we may neglect T 2 in the de- 
nominator. We then find 



Quu = p^C 2 A G 2 F n B T b T 

4 7T 

= 2.37 x 10 35 erg cm" 3 s" 1 



x Pl4 ^10 



2 + Tio ' 



(CIO) 



where we have used Eq. (C7) for the spin re- 
laxation rate, and Tio = T/10 MeV and pu = 
p/10 14 g cm~ 3 . 

C.3. Neutrino Absorption 

We next ask for the neutrino mean free path 
against the absorption process vDNN — ^ NN (in- 
verse bremsstrahlung). By means of the structure 



function one easily finds for a neutrino of energy e 



brcms 



C\G 2 1 
— nB 2e 



d 3 k 



2e(27r) 3 



f(e)24eeS(e+e), 



(Cll) 



where over-barred quantities refer to the anti- 
neutrino which is absorbed together with our test 
neutrino. Using Boltzmann statistics and assum- 
ing thermal equilibrium, the distribution function 
is /(e) = e -£ / T . Ignoring T 2 downstairs in the 
Lorentzian structure function, we find 



brcms 



3C%G%n B r 2 r 2 
ir 2 5 e 



B(e/T) 



1.92 km" 1 J—^— B(e/T). 

eio V 2 + T 10 V ' ' 



(C12) 



The dimensionless function 



B(x 



dx 



(x + x) 2 e x + e +x 



(C13) 



is shown in Fig. 23. For neutrino energies of a few 
T, corresponding to x-values of a few, this function 
is close to unity. It varies slowly for most energies 
of interest. Therefore, for rough estimates we may 
use that the bremsstrahlung absorption rate varies 
as e _1 . 




Fig. 23.— Function B{x) defined in Eq. (C13). 
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C.4. Neutrino Emissivity 

We also need the "source strength" or emissiv- 
ity of the medium in neutrinos. In contrast with 
the energy-loss rate, we need the production rate 
of neutrinos, not of pairs. Ignoring Pauli-blocking 
(we always use Boltzmann statistics!), the differ- 
ential production rate can be found directly from 
an expression like Eq. (CI) or else by detailed- 
balance arguments from the bremsstrahlung ab- 
sorption rate, leading to 

to") 

As the absorption rate scales roughly as e _1 , the 
production rate varies approximately as ee~ e ^ T . 



D. Recoil and Inelastic Scattering Com- 
bined 

The dynamical structure function for brems- 
strahlung also accounts for inelastic ncutrino- 
nucleon scattering, v + NN — > NN + v, i.e. 
neutrinos can gain or lose energy in a collision 
with a nucleon due to the simultaneous interac- 
tion of the nucleon with another nucleon. This 
energy gain or loss is unrelated to nucleon recoil. 
One could ignore nucleon recoils, and still take 
neutrino-nucleon energy transfers into account by 
using the Lorentzian structure function Eq. (C2) 
in the differential scattering cross section Eq. (Bl). 
However, while one can have recoil effects without 
the bremsstrahlung-related inelastic energy trans- 
fer if the medium is sufficiently dilute, one cannot 
have inelastic scattering without recoil as both ef- 
fects arise when the nucleon mass is not taken to 
be infinite. 

In principle, there is one complete dynamical 
structure function which includes all effects, but 
alas, it is not known. In principle, it should be pos- 
sible to calculate a dynamical structure function 
which accounts for both bremsstrahlung and recoil 
effects. However, all bremsstrahlung calculations 
were performed in the long- wavelength limit where 
the momentum transfer of the neutrinos to the nu- 
cleon pair was ignored. In this limit, recoil effects 
do not show up. 

We presently derive a heuristic form for S(u, k) 
as a convolution of the expressions for recoil and 
bremsstrahlung. The most important property 



of the structure function is detailed balancing 
Eq. (B4). As this property would not persist after 
convolving two different structure functions with 
each other, we use the symmetric form S(lu, k) as 
a starting point. The proper structure function is 
then trivially recovered by Eq. (B8). 

The symmetric form of the recoil structure 
function is explicitly 



5 rccoil (^,fc) = 4 /JL e -<- 2 /4Tu, fce - Wfc /4T 



uj k T 

w/2T _|_ e -^/2T 



(Dl) 



In the non-relativistic limit of a very large nu- 
cleon mass m we have wt < T so that the en- 
ergy transfers are typically much smaller than T. 

Therefore, in this limit we may approximately use 
e -u, fc /4T = e ±^/4T = l SQ that the stl . ucture f unc . 

tion simplifies to 



Srecoil(u') fc) — — 



1 271171 



exp 



m /wv 
~2T \k) 



(D2) 

In the limit of m — > oo the two forms are equiv- 
alent. The new version gives an average energy 
transfer per collision which is 15% smaller than 
the original version for T = 10 McV. The original 
version was derived under the assumption of non- 
relativistic nucleons even though for T = 10 MeV 
a typical nucleon speed is about 0.17 of the speed 
of light. Therefore, at the level of approximation 
where nucleons are treated as non-relativistic it is 
not even clear if the simplified version or the orig- 
inal one are a better approximation to the true 
state of affairs. We have introduced the simpli- 
fied symmetric version primarily because it allows 
us to combine it easily with inelastic neutrino- 
nucleon scattering. 

Bremsstrahlung is characterized by the Loren- 
tzian structure function Eq. (C2), which is itself a 
heuristic result. For the present purpose it is more 
useful to replace it with the ansatz 



'brcms 



-(*/4) ( w /r) 



') . (D3) 



Crucially this function has the same behavior 
2T/u> 2 for w > T and it obeys the same normal- 
ization. 

After some tinkering one can derive an interpo- 
lation formula which roughly represents the con- 
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volution of the two expressions, 



2^ 



k + 7eT/ K 



+ 27 



1 — exp < — — 



4 \-f + Ke-~'/ K 



(D4) 



which is properly normalized, even though it may 
not look that way. The r.h.s. is expressed in terms 
of the dimensionless variables 

x = — , 
T 

k Fxf 
r 

T' 



7 = 



(D5) 



The interpolation formula takes on the relevant 
limiting cases for 7 — > (no bremsstrahlung) or 
k — ► (no recoil). 

The distribution dajdti of final-state neutrino 
energies is again obtained by integrating Eq. (Bl) 
over cos# with Eq. (D4) for S(w,k) and observing 
that w = ei — e 2 and k = (ef + t\ - 2eie 2 cos#) 1//2 . 
In analogy to Fig. 2, we show in Fig. 24 this 
distribution for T = 10 MeV and an initial- 
state neutrino energy ei = 30 MeV for sev- 
eral values of the assumed spin relaxation rate 




20 40 60 80 

e 2 [MeV] 

Fig. 24. — Distribution of final-state energies e 2 
of a neutrino with initial energy ei = 30 MeV, 
scattering on non-degenerate nucleons in thermal 
equilibrium with T = 10 MeV. As a structure 
function, Stotai of Eq. (D4) was used, including 
both nucleon recoil and inelastic scattering. From 
top to bottom, the curves are for 7 = 0, 0.5, 1.0, 
1.5 and 2.0. 



7. We note that 7 = 1.0 corresponds roughly 
to a density of 10 14 g cm~ 3 . Therefore, while 
the bremsstrahlung-related broadening of the re- 
coil structure function can be quite significant, this 
is only the case for relatively large densities. 



E. Numerical Code 

In order to study and v T transport we have 
written a simple Monte Carlo code to solve the 
Boltzmann collision equation for a plane-parallel 
geometry. We prescribe a power-law density and 
temperature profile, i.e. in contrast to full-scale SN 
codes there are no "radial" zones. The medium 
consists of one species of non-degenerate nucle- 
ons with an adjustable mass. The only micro- 
scopic processes are vN scattering and NN brems- 
strahlung. We treat bremsstrahlung as described 
in Appendix C. In particular, to calculate the 
mfp against inverse bremsstrahlung absorption 
uuNN — > NN we assume that the other neutrino 
is in thermal equilibrium. In order to include en- 
ergy transfer in vN collisions we use the dynamical 
structure function described in Appendix D where 
the effect of nucleon recoil and inelastic vN scat- 
terings are combined. We can switch recoil and 
inelastic scattering independently on or off. 

To solve for the neutrino distribution function 
we follow the trajectory of a single neutrino, allow- 
ing for random collisions or absorptions. When- 
ever the neutrino is lost through the surface, into 
the blackbody surface at the bottom, or when it 
has been absorbed by inverse bremsstrahlung, a 
new neutrino is launched. Its starting point can be 
somewhere in the medium according to the source 
strength of the NN bremsstrahlung process, or it 
can be launched at the bottom from the blackbody 
surface. If NN bremsstrahlung is switched off, 
the latter is the only option. When bremsstrah- 
lung is on, we can switch the blackbody source 
off, creating all neutrinos within the medium, or 
we may replace neutrinos absorbed by the black- 
body surface with neutrinos emitted by it. When 
TV TV bremsstrahlung is on, the neutrinos achieve 
full thermal equilibrium below the energy sphere; 
the treatment of the lower boundary condition is 
then found to be irrelevant as long as it is some 
distance from the energy sphere. 

This approach makes it impractical to include 
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neutrino phase-space blocking effects as one would 
need the distribution function as input informa- 
tion for the scattering rates. Since phase-space 
blocking effects are not important for fx- and r- 
neutrinos in the atmosphere of a SN core, we have 
opted for ignoring Pauli blocking entirely and to 
treat neutrinos with Maxwcll-Boltzmann statis- 
tics. Put another way, the equilibrium distribu- 
tion found by our code is e~ f -l T , not a Fermi-Dirac 
distribution. 

As the neutrinos propagate, we sample their en- 
ergy and direction whenever they cross any of a set 
of sampling radii that were chosen at the begin- 
ning of the run. This way we accumulate statistics 
for the distribution function and several of its an- 
gular and energy moments at the sampling radii. 
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